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ABSTRACT 


The  air  velocity  profile  of  a  U.  S.  Army  Research  Institute  of  Environmental 
Medicine  (USARIEM)  test  chamber  was  mapped  using  an  anemometer  tree 
system.  This  particular  chamber  houses  an  articulated  copper  manikin.  It  was 
found  that  the  air  velocity  within  the  chamber  is  a  disparate  property.  Airflow 
patterns  are  significantly  different  at  different  chamber  fan  speed  settings,  and  the 
practice  of  employing  a  single  thermal  or  cup  anemometer  to  monitor  chamber  air 
velocity,  in  general,  is  inadequate.  It  is  recommended  that  the  air  velocity  profile 
of  a  given  test  chamber  should  be  mapped  by  techniques  described  in  this  report, 
for  all  protocol  studies  in  which  air  velocity  data  are  important,  such  as  clothing 
insulation  evaluation  and  convective  heat  transfer  studies. 


v 


INTRODUCTION 


Selection  of  proper  clothing  has  an  important  influence  on  the  performance 
of  military  personnel.  Biophysical  models  which  allow  prediction  of  clothing 
performance  under  field  conditions  can  be  a  valuable  tool  for  development  and 
validation.  However,  all  models  are  limited  by  the  accuracy  of  the  available  test 
data.  Since  environmental  chambers  provide  proper  control  for  simulation  of  field 
conditions,  they  are  often  relied  upon  to  provide  the  initial  data  for  development 
and  the  final  data  for  validation. 

In  clothing  performance  evaluation,  the  convective  heat  exchange  is  a  critical 
variable,  and  accurate  air  velocity  data  are  necessary.  Currently,  air  velocity 
profiles  within  a  test  chamber  are,  in  general,  not  available.  In  many  test 
chambers,  the  airflow  is  assumed  to  be  uniform  throughout.  A  single  thermal  or 
cup  anemometer  is  used  to  monitor  the  air  velocity  for  the  entire  chamber.  This 
present  study  reports  on  a  technique  to  validate  the  air  velocity  scenarios.  A 
fast  responding  multi-anemometer  system  (TSI  anemometer  tree,  see  Figure  l)  was 
used  to  map  the  air  velocity  profile  of  one  USARIEM  chamber,  which  houses  an 
articulated  copper  manikin. 
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METHOD 

The  test  chamber  was  cleared  of  all  articles  except  the  articulated  manikin 
and  its  supporting  frame.  The  chamber  temperature  was  set  at  27°  C.  Chamber 
fan  speeds  were  set  to  20%,  40%,  60%  and  80%  of  maximum  fan  speed.  These 
four  settings  provided  chamber  air  velocity  ranges  between  0.2  m/s  to  2.5  m/s. 
On  the  anemometer  tree  stack,  six  omni-directional  thermal  anemometers  were 
positioned,  each  at  a  different  height  level.  The  distance  between  each 
anemometer  probe  was  30  cm,  with  the  lowest  probe  30  cm  above  the  chamber 
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floor.  The  highest  probe  was,  therefore,  180  cm  above  ground.  The  30  cm  and 
180  cm  levels  corresponded  approximately  to  the  ankle  and  the  head  positions, 
respectively,  of  the  articulated  manikin. 

Figure  2  shows  the  placements  (O)  of  the  anemometer  tree  within  the 
chamber.  In  Figure  2,  the  row  designation  of  A,  B,  C,  D,  and  E  and  the 
column  number  1,  2,  3,  4,  and  5  help  to  identify  each  tree  placement  location. 
The  anemometer  tree  was  first  placed  at  the  Al  location  in  the  chamber. 
Chamber  fan  speed  was  set  to  20%  and  air  velocity  measured.  The  chamber  fan 
speed  was  then  increased  to  the  next  level  for  the  next  measurement.  After  all 
four  wind  speeds  (20%,  40%,  60%,  and  80%)  were  completed,  the  anemometer 
tree  was  then  moved  to  another  location,  to  begin  another  round  of  measurement. 
The  sequence  of  anemometer  tree  placement  was  Al,  A2,  A3,  A4,  A5,  Bl,  B2,  ... 
B5,  Cl,  ...  C5,  Dl,  ...  D5,  El,  ...  E5. 

Wind  speed  was  measured  using  a  portable  (COMPAQ)  computer.  The 
thermal  anemometer  signal  was  converted  through  a  analog-to-digital  (A/D) 
converter,  to  digital  data.  The  digital  data  were  then  sampled  and  stored  in  the 
computer.  Figure  3  shows  a  schematic  diagram  of  the  signal/data  path.  At 
each  tree  location,  for  each  chamber  fan  speed  setting,  the  anemometer  data  were 
collected  continuously  over  a  thirty  minute  period.  The  measured  air  velocity 
shown  in  Tables  1-4  represent  the  average  of  the  30  minute  period.  With 
twenty-five  tree  placements,  each  at  four  different  chamber  wind  speeds,  a  total  of 
100  sets  of  air  velocity  profile  data  were  collected. 

RESULTS 

The  entire  set  of  the  measured  anemometer  data  are  included  as  Tables  1 
through  4.  Tables  1,  2,  3  and  4  contain,  respectively,  the  air  velocity  data  at 


20%,  40%,  60%  and  80%  of  maximum  chamber  fan  speed.  The  anemometer  tree 
locations  are  designated  using  row  A,  B,  C,  D,  and  E,  and  column  1,  2,  3,  4, 
and  5  (refer  to  Figure  2).  The  air  velocity  profiles  are  plotted  in  three 
dimensional  maps,  shown  in  Figures  5  through  10.  Figure  4  is  used  as  an 
example  to  explain  the  map  schematics.  The  chamber  floor  is  represented  by  a 
rectangle.  Direction  of  airflow  is  from  right  to  left.  On  the  chamber  floor,  each 
anemometer  tree  placement  location  is  designated  by  a  ’+’  symbol.  The  human 
figure  represents  the  position  of  the  articulated  manikin.  The  manikin  is  deleted 
in  most  plots  to  provide  a  clearer  background  for  the  air  velocity  map.  The 
level  label  (120  cm)  indicates  the  height  level  of  the  anemometers.  The  air 
velocity  data  are  mapped  in  a  grid  plot.  The  grid  in  Figure  4  represents  the  air 
velocity  map  from  anemometer  measurement  at  the  60  cm  level  (at  80%  fan 
speed).  If  the  air  velocity  of  a  particular  level  is  fairly  uniform,  then  a  flat  grid 
is  the  result.  An  uneven  grid  map  indicates  large  variations  in  the  air  velocity 
pattern.  The  maximum  and  minimum  air  velocity  values  for  the  data  set 
represented  in  the  grid  are  also  displayed  to  show  the  range.  The  two  digit 
number  at  the  upper  right  hand  corner  gives  the  fan  speed  and  anemometer  level. 
The  first  digit,  range  1  to  4,  gives  the  fan  speed:  1  =  20%,  2  =  40%,  3  =  60%, 
and  4  =  80%  of  maximum  chamber  fan  speed.  The  second  digit,  range  1  to  6, 
indicates  the  anemometer  level:  1  =  30  cm,  2  =  60  cm,  ...  6  =  180  cm.  The 
anemometer  level  is,  of  course,  already  explicitly  displayed. 

Figures  5-10  display  the  full  set  of  air  velocity  data  contained  in  Tables  1 
-  4.  Each  of  the  six  figures  represent  one  anemometer  height  level.  Within  each 
figure,  four  plots  <a>,  <b>,  <c>,  and  <d>  are  derived  from  the  four  fan  speed 
settings  20%,  40%,  60%,  and  80%,  respectively. 


The  computer  algorithm  used  to  produced  the  simulated  three  dimensional 
air  velocity  map  is  included  as  Appendix.  The  program  was  written  using 
QuickBASIC  Compiler  (version  4.0).  It  can  be  easily  converted  to  run  on  a 
BASIC  interpreter. 

DISCUSSION 

Several  observations  can  be  made  from  Figures  5  -  10.  One  immediately 
apparent  is  that  the  variation  in  air  velocity  was  much  more  pronounced  at 
higher  fan  speed  than  at  lower  fan  speed.  Plots  <a>  in  Figures  5-10  show  at 
20%  fan  speed,  air  velocity  was  nearly  uniform  throughout  the  chamber,  at  all 
height  levels.  It  should  be  mentioned  that  the  20%  fan  speed  is  the  minimal  air 
circulation  required  to  maintain  chamber  temperature  regulation.  In  contrast, 
plots  <d>  of  Figures  5-10  displayed  large  variation  in  air  velocity  when  80% 
chamber  fan  speed  was  employed. 

The  second  observation  is  that,  as  expected,  structures  within  the  chamber 
alter  significantly  the  airflow  pattern.  This  includes  all  test  equipment  used  in 
the  chamber  and  other  built-in  structures  and  constructions  within  the  chamber. 
In  this  particular  testing  chamber,  the  articulated  manikin  and  its  supporting 
frame  are  the  most  prominent  structures.  In  Figure  8,  at  120  cm,  it  can  be  seen 
that  air  velocity  decreased  not  only  behind  but  also  in  front  of  the  manikin.  In 
comparison  to  Figure  10,  where  the  anemometer  height  (180  cm)  was  slightly 
above  the  head  of  the  manikin,  interference  of  the  airflow  was  much  less 
noticeable. 

There  were  also  variations  in  the  air  velocity  maps  which  are  not  as  easily 
explained.  In  Figure  9,  plot  <d>,  aside  from  the  well  around  the  manikin,  there 
were  two  other  dips  in  the  air  velocity  grid:  at  positions  B2  and  D3  (refer  to 
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Figure  2  for  position  designation).  Also,  dips  at  positions  B2  and  E2  in  Figure 
8<d>,  and  at  positions  B2  and  D4  in  Figure  7<d>,  were  evident.  These 
decreases  might  be  caused  by  the  manikin  supporting  frame,  although  the  evidence 
is  not  conclusive.  If  the  supporting  frame  were  the  culprit,  symmetrical  effect 
should  have  been  found,  i.e.,  the  dip  at  B2  should  couple  with  another  dip  at 
D2,  since  the  frame  is  symmetrical  about  the  manikin.  No  symmetry  was 
evident,  and  no  conclusion  can  be  made.  Nevertheless,  these  variations  must  be 
accepted  as  inherent  to  the  particular  chamber  construction,  and  their  effects  must 
be  properly  considered  when  using  any  given  chamber. 

In  a  recent  study  in  another  environmental  chamber,  the  USARIEM  High 
Elevation  Simulation  Facility  (altitude  chamber),  the  authors  noticed  that  not  only 
does  the  chamber  exhibit  a  characteristic  airflow  pattern  but  also  variants  in 
airflow  can  be  attributed  to  different  simulated  elevations.  For  example,  during 
the  ascending  and  descending  phases,  the  airflow  showed  marked  variation  during 
a  transient  period  before  slowly  stabilizing  to  a  new  characteristic  pattern  of  the 
new  chamber  environment. 

It  seems  clear  that  no  environmental  chamber  should  be  assumed  to  exhibit 
a  laminar  or  uniform  flow  pattern.  Each  individual  chamber  has  its  own  inherent 
air  velocity  profile,  perhaps  due  to  its  construction,  and  the  profile  invariably 
changes  as  the  chamber  condition  is  altered.  The  standardized  technique 
described  in  this  report  could  be  easily  employed  prior  to  or  along  with  specific 
studies  requiring  precise  estimation  of  heat  transfer  coefficients  attributed  to 
airflow  variations. 


CONCLUSION 


The  USARIEM  articulated  manikin  chamber  air  velocity  profile  was  mapped 
using  an  anemometer  tree  system.  It  was  found  that  the  air  velocity  within  the 
chamber  is  far  from  uniform,  and  the  airflow  patterns  can  be  significantly 
different  at  different  chamber  fan  speed  setting.  As  expected,  equipment  and 
structure  within  the  chamber  do  alter  the  airflow  pattern.  Furthermore,  there  are 
also  variations  in  the  airflow  which  are  apparently  inherent  to  the  chamber 
construction  itself,  but  can  not  otherwise  be  explained.  The  assumption  of 
uniform  airflow  within  a  chamber  is  inaccurate  and  the  practice  of  employing  a 
single  anemometer  to  monitor  the  chamber  air  velocity  is  inadequate.  It  is 
recommended  that  an  air  velocity  profile  of  given  chamber  should  be  mapped  for 
all  protocol  studies  in  which  air  velocity  and/or  precise  heat  transfer  coefficients 
are  important. 


Table  l 


180cm 


0.226 

0.274 

0.287 

0.302 

0.258 

0.325 

0.287 

0.183 

0.225 

0.201 

0.307 

0.298 

0.252 

0.237 

0.178 

0.258 

0.283 

0.268 

0.205 

0.256 

0.163 

0.172 

0.259 

0.257 

0.288 


Air  velocity  data  (m/s)  at  20%  maximum  fan  speed 
refer  to  Figure  2  for  anemometer  location  designation 


20%  raaxinnxn  fan  speed 


Anemometer  Probe  Height 


Anemome  t  e  r 


150cm  120cm  90cm  60cm 


30cm 


Location 


0 

239 

0 

264 

0 

235 

0 

295 

0 

274 

A1 

0 

264 

0 

343 

0 

314 

0 

375 

0 

327 

A2 

0 

286 

0 

349 

0 

348 

0 

346 

0 

348 

A3 

0 

246 

0 

336 

0 

325 

0 

366 

0 

317 

A4 

0 

272 

0 

322 

0 

350 

0 

336 

0 

364 

A5 

0 

254 

0 

349 

0 

339 

0 

361 

0 

334 

B1 

0 

151 

0 

238 

0 

164 

0 

287 

0 

229 

B2 

0 

126 

0 

249 

0 

205 

0 

311 

0 

264 

B3 

0 

169 

0 

184 

0 

245 

0 

295 

0 

265 

B4 

0 

194 

0 

113 

0 

187 

0 

241 

0 

260 

B5 

0 

246 

0 

302 

0 

247 

0 

288 

0 

215 

Cl 

0 

241 

0 

285 

0 

234 

0 

286 

0 

199 

C2 

0 

250 

0 

304 

0 

301 

0 

340 

0 

249 

C3 

0 

148 

0 

216 

0 

208 

0 

280 

0 

210 

C4 

0 

110 

0 

070 

0 

210 

0 

163 

0 

221 

C5 

0 

270 

0 

305 

0 

291 

0 

303 

0 

305 

D1 

0 

207 

0 

242 

0 

180 

0 

262 

0 

229 

D2 

0 

060 

0 

249 

0 

105 

0 

107 

0 

100 

D3 

0 

177 

0 

.185 

0 

095 

0 

105 

0 

148 

D4 

0 

089 

0 

165 

0 

166 

0 

166 

0 

167 

D5 

0 

169 

0 

194 

0 

230 

0 

321 

0 

357 

El 

0 

183 

0 

231 

0 

313 

0 

409 

0 

428 

E2 

0 

342 

0 

402 

0 

406 

0 

385 

0 

408 

E3 

0 

387 

0 

385 

0 

299 

0 

263 

0 

262 

E4 

0 

288 

0 

284 

0 

262 

0 

227 

0 

217 

E5 

7 


t 
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Table  2  Air  velocity  data  (m/s)  at  40%  maximum  fan  speed 

refer  to  Figure  2  for  anemometer  location  designation 


Table  3 


180cm 


1.262 

1.513 

1.556 

1.546 

1.366 

1.402 

1.452 

1.094 

0.849 

0.698 

1.530 
1.597 

1.531 
1.107 
0.852 
1.280 
1.628 
1.378 
1.275 
1.133 
0.591 
0.699 
1.034 
0.837 
0.924 


Air  velocity  data  (m/s)  at  60%  maxinun  fan  speed 
refer  to  Figure  2  for  anemometer  location  designation 


60%  maxirun  fan  speed 


Anemometer  Probe  Height 


Anemometer 


150cm  120cra  90cm  60cm 


30cm 


Local  ion 


1 

096 

1.625 

1 

245 

A1 

1 

365 

1.896 

1 

539 

A2 

1 

303 

1 

676 

1 

630 

1.789 

1 

515 

A3 

1 

324 

1 

658 

1 

621 

1.707 

1 

451 

A4 

1 

258 

1 

553 

1 

533 

1.458 

1 

408 

A5 

0 

996 

1 

547 

1 

301 

1.579 

1 

139 

B1 

0 

724 

1 

004 

0 

698 

1.264 

0 

965 

B2 

0 

558 

1 

223 

1 

049 

1.602 

1 

098 

B3 

0 

740 

0 

696 

0 

920 

1.562 

1 

001 

B4 

0 

978 

0 

489 

0 

974 

1.127 

0 

974 

B5 

1 

229 

1 

360 

1 

242 

1.760 

1 

057 

Cl 

1 

395 

1 

288 

1 

253 

1.909 

1 

185 

C2 

0 

799 

1 

379 

1 

410 

1.712 

1 

048 

C3 

0 

642 

0 

956 

1 

028 

1.451 

0 

862 

C4 

0 

379 

0 

359 

0 

911 

0.805 

0 

934 

C5 

1 

251 

1 

320 

1 

282 

1.470 

1 

281 

D1 

1 

276 

1 

288 

1 

185 

1.309 

1 

115 

D2 

0 

358 

1 

214 

0 

759 

1.304 

0 

575 

D3 

0 

574 

0 

883 

0 

630 

1.029 

0 

863 

D4 

0 

532 

0 

989 

0 

884 

1.082 

0 

942 

D5 

0 

764 

0 

659 

0 

614 

0.960 

0 

811 

El 

0 

721 

0 

595 

0 

531 

0.740 

0 

718 

E2 

1 

210 

1 

080 

0 

853 

1 .017 

1 

140 

E3 

1 

536 

1 

444 

1 

132 

1 . 197 

1 

297 

E4 

1 

599 

1 

514 

1 

233 

1  . 162 

1 

174 

E5 
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Table  4  Air  velocity  data  (m/s)  at  80%  maximun  fan  speed 

refer  to  Figure  2  for  anemometer  location  designation 


80%  maximum  fan  speed 


Anemometer  Probe  Height 


Anemome  ter 
Locat ion 


Figure  1  Anemometer  tree  (modified  from  TSI  contract  manual) 


Figure  2  Articulated  mannikin  chamber  schematic  diagram 
with  the  mannikin  and  supporting  frame 
anemometer  tree  placement  location  (#) 
direction  of  air  flow  is  from  right  to  left 


air 

flow 

direction 
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Figure  3  Diagram  showing  signal/data  path  from  the  anemometer  tree  to  the  computer 


APPENDIX 


This  program  plots  a  simulated  three  dimensional  air  velocity 
map  from  data  gathered  using  the  Anemometer  Tree  System  In  an 
environmental  chamber.  The  entire  data  set  collected  from  the 
articulated  manikin  chamber  are  included  in  this  report  as 
Tables  1-4,  and  are  deletsd  from  this  program  to  save  space. 

A  few  examples  of  the  data  set  are  left  to  show  the  Input  data 
format.  This  program  is  written  using  QuickBASIC  Compiler 
(•version  4.0).  It  can  be  converted  to  run  on  BASIC  Interpreter 
by  adding  line  numbers,  and  replacing  all  the  Labels  after  a 
GOTO  statement  with  the  appropriate  line  numbers. 


'  chamber  level  labels 

DATA  "180cm" . "150cm" , "120cm" , "90cm" , "60cm" , "30cm" 


DATA  1 .945. 1.592, 1 .845, 1 .590,2.394, 1 .766 
•  DATA 
'DATA 
'  DATA 


air  velocity  data 
(80%  fan  speed) 
see  data  set  in  Table  4 


DATA  1 . 262 , 1 . 096 . 1 . 255 . 1 . 063 . 1 . 625 , 1 . 245 
'DATA 
'  DATA 
'DATA 


air  velocity  data 
(60%  fan  speed) 
see  data  set  in  Table  3 


DATA  0.609,0.579,0.670,0.550.0.818.0.671 

'DATA 

'DATA 

'DATA 


air  velocity  data 
(40%  fan  speed) 
see  data  set  in  Table  2 


DATA  0.226,0.239,0.264,0.235,0.295,0.274 
'  DATA 
'  DATA 
'  DATA 


air  velocity  data 
(20%  fan  speed) 
see  data  set  In  Table  1 


DIM 

DIM 

DIM 

DIM 

DIM 

DIM 


Manp( 14944)  AS  INTEGER  ' 
Bord( 14102)  AS  INTEGER  ' 
Cham(4,  5.  5,  6)  AS  SINGLE 
I  eve  I  $06)  ' 

wmax(4,  6),  wmin(4,  6)  ' 


mannikin  plot  storage  array 
screen  border  plot  storage  array 
air  velocity  data  array 
anemometer  level  labels 
max  &  min  air  velocity 


vx(5,  5),  vy(5,  5)  AS  INTEGER 


air  velocity  map  points 


SCREEN  9 
CLS 


backgd  -  1  '  mannikin  background  Indicator  (1-display,  0-no) 


px  -  375:  py  -  305  ' 
ax  -  25:  ay  -  80  ' 
lx  -  270:  ly  -  270  ' 
SX  -  50:  sy  -  55  ' 
dx  -  30:  dy  -  75  ' 


pset  (starting)  location 
turn  angle  DRAW  degrees 
border  DRAW  statement  values 
probe  location  DRAW  values 
wind  direction  arrow  angle  DRAW  degrees 


14 


draw  mannikin 


CIRCLE  (176,  85),  14,  , 
DRAW  “AO  8015”  +  “TA«“  + 
DRAW  "TA-"  +  VARPTRS ( ax ) 
DRAW  "TA-“  ♦  VARPTRS (ax) 
DRAW  “TA-”  +  VARPTRS(ax) 
DRAW  "TA-"  ♦  VARPTR$(ax) 
DRAW  "TA-"  +  VARPTR$ ( ax ) 


.  16  /  15  ' 

VARPTR$(ax)  +  "R3 1  AO  D70 


+ 

"LI  1 

AO 

NU30  D70 

♦ 

"L20 

AO 

NU60" 

+ 

“L20 

AO 

U70  NU30 

♦ 

"LI  1 

AO 

U70" 

♦ 

“R31 

M 

GET  (70,  70)-( 560,  310),  Manp 


SAVE  mann Ik i n 


CLS 

PSET  (px,  py)  *  draw  probe  locations 

FOR  y  -  1  TO  5 

FOR  x  -  1  TO  5 

DRAW  “AO  NU2  ND2  NR2  NL2 " 

DRAW  "TA-"  +  VARPTRS ( ax )  +  "BR«“  ♦  VARPTRS OSX) 

NEXT  x 

PSET  (px,  py):  yy  •  sy  *  y 

DRAW  "TA-"  +  VARPTRS(ay)  +  “BU-"  ♦  VARPTR$ ( yy ) 

NEXT  y 
« 


• 

draw  border  outline 

PSET  (px 

♦  4.  py 

♦  15)  ' 

draw  front  border 

FOR  I 

-  1  TO 

8  ' 

thick  x  side 

DRAW  "TA-" 

+  VARPTRSOax) 

+  "NR-" 

+  VARPTRS (lx)  +  "AO  BD1 " 

NEXT 

1 

PSET  (px 

+  4,  py 

♦  15)  ' 

right  side 

DRAW 

"TA-"  + 

VARPTRS ( ax ) 

♦ 

"8R-"  + 

VARPTRS ( lx) 

DRAW 

"TA-"  + 

VARPTRS ( ay ) 

+ 

"U270" 

PSET  (px 

♦  4.  py 

+  15) 

FOR  J 

-  1  TO 

8  ' 

thick  y  side 

DRAW  "TA-" 

+  VARPTRS  (>ay) 

+  "NU-" 

+  VARPTRS ( 1 y )  +  “AO  BD1 " 

NEXT 

J 

PSET  (px 

+  4,  py 

+  15)  * 

back  side 

DRAW 

"TA-"  + 

VARPTRS(ay) 

♦ 

"BU-"  ♦ 

VARPTRS ( ly) 

DRAW 

t 

"TA-"  + 

VARPTRS (ax) 

+ 

"R270" 

dy3  -  dy 

-  3  ' 

draw  wind  direction  arrow 

PSET  (px 

♦  192,  py  +  6) 

DRAW  "TA-"  +  VARPTRS(dy)  +  "U18"  +  "TA-"  +  VARPTR$(dx)  +  “L10“ 
PSET  (px  +  192,  py  +  6) 

DRAW  "TA-"  +  VARPTRS(dx)  +  "R40"  +  "TA-”  +  VARPTR$(dy3)  +  “U21“ 
DRAW  "TA-"  +  VARPTRS(dx)  +  "R12  TA90  U37  TA180  U26” 

LOCATE  24,  63:  PRINT  "wind  direction"; 

f 

GET  (0,  180)-(639,  349),  Bord  '  save  border  outline  &  probe  loc. 

« 

PUT  (70,  70),  Manp,  OR  '  put  back  mannikin 

FOR  z  -  6  TO  1  STEP  -1 

READ  level$(z)  '  read  In  anemometer  level  labels 

NEXT  z 
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Figure  4  Illustration  of  air  velocity  map 
( 120cm  height  level ) 


6 


wind  speed  -  4;  chamber  dim. 
w 

highest  wind  setting  -  4 
lowest  wind  setting  -  1 


5*5;  anemometer  level  - 
x  *  y  z 

6  -  180  cm 
1  -  30  cm 


FOR  w  «  4  TO  1  STEP  -1  *  for  4  chamber  wind  speed  settings 

FOR  Z  -  6  TO  1  STEP  -1 
wmax(w,  z)  -  0 

wmin(w,  z)  -  1000  *  Initialize  max  &  min  air  velocity 


NEXT  Z 


read  in  chamber  air  velocity  data 
FOR  x  -  1  TO  5  '  chamCw,  x,  y,  z) 

FOR  y  -  1  TO  5  '  chamber  dim.  »  x  *  y 

FOR  z  -  6  TO  1  STEP  -1  ’  anemometer  level  -  z 

READ  cham(w.  x,  y,  z)  ’  chamber  wind  speed  -  w 
IF  cham(w,  x,  y,  z)  >  wmax(»w,  z)  THEN 
wmax(w.  z)  -  cham(w,  x,  y,  z) 


END  IF 


IF  cham(w,  x,  y,  z)  <  wmln(iw,  z)  THEN 
wmln(w,  z)  -  cham(w,  x,  y,  z) 

END  IF 


NEXT  z 
NEXT  y 
NEXT  X 
NEXT  w 


ScrLevel :  '  ml  sc  input  &  screen  displays 

LOCATE  25.  Is  PRINT  SPC(79); 

LOCATE  25.  1 

PRINT  “Enter  wind  speed  (1-4),  chamber  level  (1-6), 

INPUT  “(11  to  46),  0  to  quit  “;  ppX 


SELECT  CASE  ppX; 

CASE  0  1  Quit 

SCREEN  0:  CLS  :  END 
CASE  11  TO  16 

wx  -  1  '  wX  -  selected  chamber  wind  speed 

z%  -  ppX  -  wX  *  10  '  z%;  -  selected  chamber  level 

CASE  21  TO  26 
wX  -  2 

z%  -  ppX  -  wX  *  10 
CASE  31  TO  36 

wX  »  3 

zX  -  ppX  -  wX  *  10 
CASE  41  TO  46 
wX  -  4 

ZX  -  ppX  -  wX  *  10 
CASE  ELSE 
BEEP 

LOCATE  24,  1;  PRINT  SPC(79); 

LOCATE  24 ,  1 

PRINT  "acceptable  Input:  11-16,  21-26,  31-36,  41-46"; 
GOTO  ScrLevel 
END  SELECT 
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— * 

1 


CLS  '  start  displaying  air  velocity  map 

PUT  (0,  180),  Bord  '  display  border  outline 

IF  backgd  «  1  THEN  PUT  (70,  70),  Manp  '  display  mannikin 

LOCATE  20  -  (2  *  ZX  +  3 ) ,  1:  PRINT  “level  levelS(zX) 

LOCATE  1,  2 

PRINT  "max  air  - 


PRINT 

USING  “*.##•"; 

wmax(wX,  zX) ; 

:  PRINT  "  m/s“ 

LOCATE  2. 

2 

PRINT 

“min  air  ■  “ ; 

PRINT 

USING  “#.#«*“; 

wmln(wX,  zX) ; 

:  PRINT  "  m/s" 

LOCATE  1 , 

• 

77:  PRINT  ppXi; 

i 

print  Input  selection 

PSET  (px. 

py)  * 

Compute  air  velocity  map 

range  -  (wmax(w%,  zX)  -  wmln(w%i,  zX)) 

FOR  y  -  1  TO  5 
FOR  x  -  1  TO  5 

tt  ■  (cham(wX,  5  -  x  +  1.  y,  zX)  /  wmax(<wX,  zX))  *  range  *  37 
tt  -  C I  NT ( 30  *  (zX  -  1)  +  tt) 
vxOx,  y)  -  POINT(O) 
vy(x.  y)  -  POINT(I)  -  tt 

DRAW  "TA-“  +  VARPTRS(ax)  +  "BR-"  +  VARPTRS(SX) 

NEXT  X 

PSET  (px,  py):  yy  -  sy  *  y 

DRAW  "TA-“  +  VARPTRJ(ay)  +  “BU-"  +  VARPTRS ( yy ) 

NEXT  y 

FOR  y  «  1  TO  5  1  Plot  air  velocity  map  ('Simulated  3D) 

FOR  x  -  1  TO  5 
IF  x  <  5  THEN 


LINE 

(vx(x. 

y) . 

vy(x. 

y)) 

-('VX(X  +  1  , 

y). 

vy(x  ♦  1, 

y) ) 

END 

IF 

IF  y  <  5 

THEN 

LINE 

(VX(X, 

y) . 

vy(x. 

y ) ) 

—  ('Vx  ( x  ,  y  + 

i). 

vy(x,  y  + 

1 )) 

END  IF 
NEXT  X 
NEXT  y 

LOCATE  25,  1 

PRINT  " <ENTER>»new  w I nd/ I  eve  I , " ; 

PRINT  "<SPACE>-remove/add  mannikin,  <ESC>-qult"; 


i 

i 

i 

! 

I 

I 

1 

! 

> 


Nextse l 1 : 
bj  «  INKEYS 
IF  b$  -  ■"  THEN 

GOTO  NextSe I  1 

ELSE  IF  b$  -  CHRS ( 27 )  THEN  '  <ESC>  to  quit 

SCREEN  0:  CLS  :  END 

ELSE  IF  b$  -  CHRS (32)  THEN  *  <SPACE>  to  remove/add  mannikin 
PUT  (70,  70),  Manp 

IF  backgd  >  1  THEN  backgd  -  0  ELSE  backgd  •  1  '  toggle 
GOTO  NextSell  '  background  display  indicator 

ELSE  IF  bS  -  CHRS (13)  THEN  *  go  to  display 

GOTO  ScrLevel  '  new  wind  speed/chamber  level 


r-i 
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ELSE  BEEP 

GOTO  NextSel 1 

END  IF 
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